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Abstract. Non-equilibrium aging dynamics in 3D Ising spin glass Cuo.sCoo.sC^-FeCls GBIC has been 
studied by zero-field cooled (ZFC) magnetization and low frequency AC magnetic susceptibility (/ = 0.05 
Hz), where T g = 3.92±0.11 K. The time dependence of the relaxation rate S(t) = (l/H)dMzFC /dint for 
the ZFC magnetization after the ZFC aging protocol, shows a peak at a characteristic time t cr near a wait 
time t w (aging behavior), corresponding to a crossover from quasi equilibrium dynamics to non-equilibrium. 
The time t cr strongly depends on t w , temperature (T), magnetic field (H), and the temperature shift 
(AT). The rejuvenation effect is observed in both \ an d x" under the T-shift and _ff-shift procedures. The 
memory of the specific spin configurations imprinted during the ZFC aging protocol can be recalled when 
the system is re-heated at a constant heating rate. The aging, rejuvenation, and memory effects observed 
in the present system are discussed in terms of the scaling concepts derived from numerical studies on 3D 
Edwards- Anderson spin glass model. 



PACS. 75.50.Lk Spin glasses and other random magnets 
- 75.40.Gb Dynamic properties 



75.10.Nr Spin-glass and other random models 



1 Introduction 

In recent years, non-equilibrium dynamics, in particular, 
aging dynamics, of spin glass (SG) systems has been exten- 



sivelv studied theoretically [T||5||S||1||S||S|[7| 


[8"||3] and experi- 


mentally 10.11,12,13, 14, 15.16.17.18, 19 




l25U2(ill27ll28j. When the SG system is Quenched from a 



high temperature above the SG transition temperature T g 
to a low temperature T below T g (this process is called the 
zero- field cooled (ZFC) aging protocol), the initial state is 
not thermodynamically stable and relaxes to more sta- 
ble state. The aging behaviors depend strongly on their 
thermal history within the SG phase. A typical experi- 
mental method for the study of the aging dynamics is the 
time dependence of the ZFC magnetization (Mzfc) after 
the ZFC aging protocol. When a small magnetic field (H) 
is applied after an isothermal aging at a temperature T 
below T g for a wait time t w , Mzfc increases as an ob- 
served time (t) increases. The rejuvenation (chaos) and 
memory effects are also significant features of the aging 
dynamics. These effects are typically measured from the 
low frequency AC magnetic susceptibility. 

Such an aging dynamics of SG phase is explained mainly 
in terms of a real-space picture (the droplet model) pQ. 
In this picture, the SG coherence length for equilibrium 
SG order, grows up slowly in aging processes. The scal- 
ing properties of age-dependent macroscopic susceptibil- 
ity can be described by a growing coherence length Lx(t). 



The droplet model also predicts the following two. (i) The 
equilibrium SG states at two temperatures with the dif- 
ference AT are uncorrelated in length scales larger than 
the overlap length L^t, i.e., so-called temperature (T)- 
chaos nature of the SG phase, (ii) In the equilibrium and 
thermodynamic limits the SG phase is broken by a static 
magnetic field H of infinitesimal strength, thereby intro- 
duced is the crossover length Lh- This length separates 
the mean size of SG domains Lx(t), such that they are 
dominated by the Zeemann energy for Lx{t) > Lh and by 
the SG free energy gap for Lx(t) <C Lh- Recently Monte 
Carlo (MC) simulations on the T- and ii-shift aging 
processes have been carried out for the three-dimensional 
(3D) Ising Edwards- Anderson (EA) SG model with Gaus- 
sian nearest-neighbor interactions with zero mean and vari- 
ance J, where J is in the unit of energy. In the T-shift 
process, only a precursor of the temperature-chaos effect is 
observed, while the results on the iTshift process strongly 
support the droplet picture that the SG phase under a fi- 
nite H is unstable in the equilibrium and thermodynamic 
limits. 

In this paper we report our results on the aging dynam- 
ics (aging, rejuvenation and memory) of a 3D Ising SG, 
Cuo.5Coo.5Cl2-FeCl3 graphite bi-intercalation compound 
(GBIC). The aging behavior of the ZFC magnetization 
has been measured at various aging processes including 
the T- and iTshift perturbations. The rejuvenation effect 
of the low frequency AC magnetic susceptibility has been 
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measured at aging processes including the T- and ii-shift 
perturbations. The memory of the specific spin configu- 
rations imprinted during the ZFC aging protocol can be 
recalled when the system is re-heated at a constant heat- 
ing rate. The existence of the characteristic lengths Lat 
and Lh is examined by the T- and i7-shift aging protocols. 
The relaxation rate defined by S(t) = (l/H)dMzFc/dixit 
(see Sec. l2.2l for more definitions) exhibits a peak at a char- 
acteristic time t cr , which is comparable to a wait time t w . 
We show that the ratio t cr /t w strongly depends on the 
aging process. These results are compared with those re- 
ported for typical spin glass systems such as Cu (1.5-13.5 

[UniEll, CdCri.7Ino.3S4 [nUH Feo.glvWsTiOa EE1 
EOj, as well as results from MC simulations pEITFIIS] . 

The equilibrium dynamics of Cuo.sCoo.sC^-FeCla GBIC 
has been reported in a previous paper . The aging dy- 
namics has been also studied: the wi-scaling of \' and 
x" is also confirmed, where uj is the angular frequency of 
the AC magnetic field. This compound undergoes a SG 
transition at T g (= 3.92 ± 0.11 K). The system shows a 
dynamic behavior that has some similarities and some sig- 
nificant differences compared to a 3D Ising SG. It shows 
critical slowing down with a value of the dynamic criti- 
cal exponent that is rather similar to an Ising SG: z = 
6.6 ± 1.2, = 0.13 ± 0.02, and \j) = 0.24 ± 0.02, where 
z, ip, and are the dynamic critical exponent, the stiff- 
ness exponent (energy exponent) which determines the 
free energy of a droplet excitation, and the barrier ex- 
ponent which describes how the barrier heights change 
with length scale. The critical relaxation time r is well 
described by r — t*(1 — T/T g )~ x with the dynamic criti- 
cal exponent x = 10.3 ± 0.7 and r* = (5.29 ± 0.07) x 10~ 6 
sec. The in-field dynamics indicates, as for an Ising SG, 
that the SG transition is destroyed by a magnetic field. 
The equilibrium dynamics shows a frequency dependence 
that is different from an Ising SG, the absorption decreases 
with increasing frequency, whereas ordinary 3D Ising SG 
shows a increasing absorption with increasing frequency. 
An aging behavior is observed that is rejuvenated by a 
large enough (magnetic field) perturbation. 

2 Background 

2.1 Scaling Properties 

We present a simple review on the aging behavior of SG 
phase after the ZFC aging protocol, based on the droplet 
model [Tir2inii] . This ZFC aging protocol process to the 
SG phase is completed at t a = 0, where t a is defined as an 
age (the total time after the ZFC aging protocol process). 
Then the system is aged at T under H = until t a = 
t w , where t w is a wait time. Correspondingly, the size of 
domain defined by Rrita) grows with the age of t a and 
reaches RT(t w ) just before the field is turned on at t = 
or t a — t w . The aging behavior in Mzfc is observed as 
a function of the observation time t. After t = 0, a prob- 
ing length Lxit) corresponding to the maximum size of 
excitation grows with t, in a similar way as R,T{t a ). When 



Lr(t) <C Rr(t w ), quasi-equilibrium dynamics is probed, 
but when Lx(t) 3> RT(t w ), non-equilibrium dynamics is 
probed. It is theoretically predicted that the mean SG 
domain-size Lt(T) is described by a power law given by 

12 

L T (<)/£o«(*Ao) lMT) , (1) 

or by P 

L T (t)/L *[(T/A g )ln(t/t )]W, (2) 

where Lq and to are microscopic length and time scale, the 
exponent l/z(T) (= bT/T g with b = 0.16) linearly de- 
pends on T except for the region near T g , and A g and ijj are 
the characteristic scale of energy barrier of droplet excita- 
tions and the associated exponent, respectively. It is pre- 
dicted by Komori et al. [2] that Ax" obeys the wt-scaling 
form given by A X " « M)" b " where b" = (d- 0)/z(T), d 
(— 3) is the dimension, and 9 = 0.20 ± 0.03. Experimen- 
tally, as shown in our previous paper |26|. Ax"(to,t) at 
3.75 K obeys the ensealing law, Ax"{u>, t) k, (uit)~ b with 
b" = 0.255 ± 0.005, where Ax"(oJ, t) = t) - (xd'M), 
where (x '(w)) is the stationary frequency-dependent ab- 
sorption. 

2.2 S{t) and x"(t w ;t + t w ) 

Here we present a simple review on the relation between 
S(t) and x" ■ The absorption x" is evaluated from the spin 
auto-correlation function C(t a — t; t a ) = {Si(t a — t)Si(t a )} 
using the fluctuation-dissipation theorem (FDT) as 0E| 

X "{At u ;t+At u ) « (-Tr/2T)dC(At u ;t+At u )/dlnt, (3) 

where t a = t + At w , At w = 2tt/uj (typically At^ < 10 2 
sec) and t is much larger than At^. In the auto-correlation 
function, the over-line denotes the average over sites and 
over different realizations of bond disorder, and the bracket 
the average over thermal noises. For slow processes, the 
dispersion x'(At u ; t + At u ) is approximated by 

X '(At u ;t + At u )*[l- C{At u ;t + AtJ)]/T. (4) 

In the quasi-equilibrium regime where the FDT holds, the 
ZFC susceptibility XzFc{tw', t+t w ) is described by 

)}/T, (5) 

where t a = t + t w and t w is a wait time: typically, t w w 
10 3 — 10 5 sec. Then the relaxation rate S(t) is described 

by 

S(t) = dxzFc(t w ;t + t w )dlnt 

= (-l/T)dC(t w ;t + t w )/dhxt, (6) 

which corresponds to (2/tt)x" (t w ; t + t w ) jl()| . 

It is predicted that C(t w ;t + t w ) can be decomposed 
into a stationary part C s t {t) and an aging part C ag it w \ t -\- 
t w ) [Sj. The latter is approximately described by a scaling 
function of Lx(t) and RT{t w ). The corresponding aging 
part of S(t) exhibits a peak at a characteristic time t cr 
{~t w ) [6Tl2"9"j . where I/r(i) ~ Rr{tw), showing a crossover 
between quasi-equilibrium region and non-equilibrium re- 
gion. 
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3 Experimental procedure 

The DC magnetization and AC susceptibility of Cuo.5Coo.5Cl2 
FeCl3 GBIC were measured using a SQUID magnetome- 
ter (Quantum Design, MPMS XL-5) with an ultra low 
field capability option. The remnant magnetic field was re- c 
duced to zero field (exactly less than 3 mOe) at 298 K for ; 
both DC magnetization and AC susceptibility measure- ■? 
ments. The AC magnetic field used in the present experi- £ 
ment has a peak magnitude of h = 0.1 Oe and frequency i 
/ = u)/2n = 0.05 Hz. Each experimental procedure for 
measurements is presented in the text and figure captions. 
The detail of sample characterizations and sample prepa- 
ration is given in the previous paper |26|. 

In our measurement of the time (i) dependence of the 
zero-field cooled (ZFC) magnetization (Mzfc), the time 
required for the ZFC aging protocol and subsequent wait 
time was precisely controlled. Typically it takes 240 ±3 sec 
to cool the system from 10 K to 3.75 K. It takes another 
t W Q — 230±3 sec until T (= 3.75 K) becomes stable within 
the uncertainty (±0.01 K). The system is kept at T = 3.75 
K and H = for a wait time t w (0 < t < 3 x 10 4 sec). 
The time for setting up a magnetic field from H = to j 
H = 5 Oe is 68 ± 2 sec. In the ZFC measurement, the 
sample is slowly moved through the pick-up coils over the 
scan length (4 cm) . The magnetic moment of the sample 
induces a magnetic flux change in the pick-up coils. It 
takes 12 sec for each scan. The data at t is regarded as 
the average of Mzfc measured over the scanning time t s 
between the times t — (t s /2) and t — (i s /2). Thus the time 
window At is a scanning time (t s ). The measurement was 
carried out at every interval of t s + t p , where t p is a pause 
between consecutive measurements. Typically we used (i) 
the time window At = 36 sec for three scans and t p = 45 
sec or 30 sec for t w (> 2.0 x 10 3 sec), and (ii) the time 
window At = 12 sec for one scan and t p = 1 or 2 sec for 
either t w (< 10 3 sec) or H > 100 Oe. 
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4 Result 

4.1 S(t) : t w and T dependence 

The t dependence of Mzfc was measured under various 
conditions. Figure H^a) shows the T dependence of the 
ZFC and FC magnetization {Mzfc and Mpc) at H = 

5 Oe, where H is applied along a direction perpendicular 
to the c axis. The system was quenched from 50 to 1.9 K 
in the absence of H before the measurement. The change 
of T with the time t during the measurement is shown 
in Fig. Bib). The ZFC magnetization was measured with 
increasing T from 1.9 to 3.75 K in the presence of H (= 
5 Oe). The system was kept at T = 3.75 K for 2.27 x 10 4 
sec. Subsequently the ZFC measurement was continued 
from 3.75 to 8 K. The system was annealed at 50 K for 

1.2 x 10 3 sec. The FC magnetization was measured from 8 
to 1.9 K. There is a remarkable increase of Mzfc during 
the one stop at T = 3.75 K. When the warming up of 
the system was restarted, the increase in Mzfc with T is 
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Fig. 1. (a) T dependence of Mzfc and Mfc at H = 5 Oe 
for Cuo.5Coo.5Ci2-FeCl3 GBIC. The measurement was carried 
out after the ZFC aging protocol: annealing of the system at 
50 K for 1.2 x 10 3 sec at H = and quenching from 50 to 
1.9 K. During the ZFC measurement the system was aging at 
3.75 K for 2.27 x 10 4 sec. (b) The change of T with t during 
the measurement of Mzfc and Mpc- The discontinuity of T 
with t at T = 4.4 K is due to the system error occurring in 
the process of temperature stabilization, (c) The relaxation 
rate S(t) (= (l/H)AM Z Fc{t)/A\rxt) at 3.75 K. The time taken 
during the measurement of Mzfc from 1.9 to 3.75 K at H — 5 
Oe was to — 5.4 x 10 3 sec. t = is a time just after T becomes 
3.75 K. 
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I l.acCj Fig. 3. T dependence of t cr at which S(t) has a peak value 

Fig. 2. (a) and (b) t dependence of S for 3.0 < T < 4.7 K. s ™,ax- t w = 2.0 x 10 3 sec. H = 1 Oe. (b) T dependence of the 

H = 1 Oe. The measurement of \zfc vs t was carried out peak height Smax (see Fig. 
after the ZFC aging protocol (annealing of the system at 50 K 
and H = for 1.2 x 10 3 sec, quenching from 50 K to T, and 
isothermal aging for t w = 2.0 x 10 3 sec), t — is the time just 

after H = 1 Oe is applied at T. time just after H = 1 Oe is applied at T. We find that 

Mzfc increases with increasing t, depending on T. Fig- 
ures |3{a) and (b) show the t dependence of S at various 
the stop: M ZFC merges with M%$ c well above T 9 . Simi- T (3.0 < T < 4.7 K), where H = 1 Oe and t u , = 2.0 x 10 3 
lar result of Mzfc vs T was reported for a 3D Ising SG sec. Each S(t) curve exhibits a peak at a characteristic 
Feo.5Mno.5Ti03. Bernardi et al. [S] have concluded that time t cr , shifting to the short t-side with increasing T. 
Mzfc is described by a scaling function of Rritw) and The peak width becomes much broader with decreasing 
L T {t) [M ZF c(T,t) = G(Lt (t) , Rt (tw ) )] j where G is the T. Figures EI a) and (b) show the T dependence of t cr 
scaling function. The t dependence of Mzfc was mon- and the peak height S max , respectively, which are derived 
itored during the stop at T — 3.75 K, where the time from Figs. EJa) and (b). The characteristic time t cr in- 
taken during the measurement from 1.9 to 3.75 K was creases with decreasing T. The T dependence of ln(t cr /t w ) 
5.4 x 10 3 sec, corresponding to a wait time t w . Figure^c) will be discussed in Sec. 0m association with the scaling 
shows the t dependence of the relaxation rate S(t), where relations. The peak height S ma x exhibits a broad peak be- 
the origin of t (t = 0) is a time when T reaches 3.75 K. tween 3.6 and 3.8 K just below T g . Similar behavior of S 
We have measured the t dependence of xzfc (= MzFc/Hfy t at various T (< T g ) has been reported by Zotev et 
at various fixed T, where t w = 2.0 x 10 3 sec and H = 1 al. for Cu (1.5 at % Mn) where S(t) is defined from 
Oe. The measurement was carried out after the ZFC aging the thermoremnant magnetization Mtr instead of Mzfc- 
protocol: annealing of the system at T — 50 K and H = The relaxation rate S(t) exhibits a peak at t — t cr , whose 
for 1.2 x 10 3 sec, quenching from 50 K to T at H — 0, and peak width becomes much broader with decreasing T. The 
isothermal aging at T for t w . The origin of t (t = 0) is a peak height S rnax decreases with decreasing T below T g . 
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Fig. 4. (a) t dependence of S for < t w < 1.5 x 10 4 sec. 
T = 3.75 K. H = 5 Oe. The ZFC aging protocol: annealing of 
the system at 50 K for 1.2 x 10 3 sec at H = 0, quenching from 
50 to 3.75 K, and then isothermal aging at 3.75 K and H = 
for a wait time t w . The measurement was started at t = 
when the field H is turned on. (b) t cr vs t TO . The straight line 
denotes a relation given by t cr = (0.68 ± 0.03)t w . (c) Scaling 
plot of S/Smax vs for the limited t w and i (5.0 x 10 < 
t w < 1.5 x 10 4 sec, < t < 6.0 x 10 4 sec). 



Figure Ufa) shows the t dependence of S at various 
t w , where T — 3.75 K and H = 5 Oe. The relaxation 
rate S(t) shows a peak at t = t cr , shifting to the long- 
t side with increasing t w . Similar behavior of S(t) vs t 
at various t w has been observed by Jonsson et al. for 
Ag (11 at.% Mn) 21 , where t cr « t w . In Fig. E^b) we 
show the characteristic time t cr as a function of t w for 
< t w < 3.0 x 10 4 sec, where a straight line denotes the re- 
lation described by t cr = (0.68±0.03)t w . Figure0Jc) shows 
the plot of S(t)/S max as a function oit/t w , where only the 
data of S vs t with long t w (5 x 10 3 < t w < 1.5 x 10 4 
sec) for < t < 6 x 10 4 sec are used. It seems that 
S{t)/ S max is well described by a scaling function of t/t w 
[S(t)/S max — F(t/t w )\ in the region of long t w , although 
the data at t w = 1.5 x 10 4 sec slightly deviates from the 
other data. The scaling function F(x) has a very broad 
peak centered at x ~ 0.68. This result suggests that the 
spin auto-correlation function C(t w , t + t w ) is described by 
a scaling function of only t/t w , since it is related to S(t) 
byEq.p. 



4.2 S(t) for < t w < 750 sec 

We have measured the t dependence of xzfc at T = 3.75 
K and H = 5 Oe as a wait time t w (0 < t w < 750 sec) is 
varied as a parameter. The measurement was carried out 
after the ZFC aging protocol: annealing of the system at 
50 K for 1.2 x 10 3 sec at H = 0, quenching from 50 to 3.75 
K, and then isothermal aging at T = 3.75 K and H = 
for a wait time t w . The origin of t (t = 0) is a time just 
after the field H is turned on. As is described in Sec. |31 
it takes t W Q — 230 ± 3 sec until the temperature becomes 
stable at 3.75 K within the experimental uncertainty of 
±0.01 K, after quenching the system from 50 to 3.75 K at 
H = 0. This time t^o is not included in the wait time t w . 
If t^o is included in t w , however, the effective wait time 
may be longer than t e J^ (= t W Q + t w ). If this is the case, 
the measurement with t w = is not possible in a strict 
sense. In spite of such a situation, here we assume that 
t is not included in t w . Figures 0Ja) and (b) show the t 
dependence of S for < t w < 750 sec, where T = 3.75 
K and H — 5 Oe. For t w = 0, S(t) shows a shoulder 
around t = 10 3 sec and a broad peak at t cr = 6.45 x 10 3 
sec. For t w — 100 sec, a small peak of S(t) is observed 
at t cr — 1.06 x 10 3 sec in addition to a possible broad 
peak around t cr = 1.0 x 10 4 sec. This small peak shifts 
to the short i-side with increasing t w : t cr = 740 sec for 
t w — 200 sec. For t w = 750 sec, a broad peak is observed 
at t cr = 2.8 x 10 3 sec. The values of t cr thus obtained are 
also plotted as a function of t w in Fig.^Jh). It should be 
noted that that for very short t w (0 < t < 200 sec), there 
are at least two kinds of domains: domains with large size 
corresponding to long t cr (« 6.5 x 10 3 sec) coexist with 
domains with small size corresponding to short t cr (« 750 
sec) in the regular aging regime. Figure shows the t 
dependence of S(t) at various T (3.4 < T < 4.5 K) and 
H = 1 Oe for t w = 0. At T = 3.2 K, S(t) shows a shoulder 
around t = 400 sec and a very broad peak at t cr = 7.1 X 10 3 
sec. This shoulder shifts to the long t-side with increasing 
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Fig. 5. (a) and (b) t dependence of 5 for < t w < 750 sec. 
T = 3.75 K. H = 5 Oe. 

T and changes in to a broad peak. At T — 3.6 K, S(t) 
shows a peak at t cr — 710 sec. It tends to increase with 
further increasing t, suggesting that t cr is longer than 10 
sec. At T = 3.8 K, S(t) shows a very broad peak centered 
at t cr = 3.8 x 10 3 sec. 

Similar behavior has been reported by Rodriguez et al. 
|25| in the time decay of the thermal remnant magnetiza- 
tion (TRM) of Cuo.94Mno.06 (T g = 31.5 K) with various 
wait time (t w = - 10 4 sec) and a series of rapid FC cool- 
ing protocol from 35 to 26 K at H (= 20 Oe). The relax- 
ation rate S(t)[= -(l/H)dM TRM [t)/dlni\ for M TRM (t) 
is equivalent to S(t) [= -(l/H)dM ZFC (t)/d\nt} for M ZFC 
They have shown that S(t) at t w — exhibits a broad 
peak at an effective time t e J* (= 19 - 406 sec), which is 
strongly dependent on the FC cooling protocol. For the 
larger t R Jf, there is a significant contamination in S(t) 
at t w (7^ 0) from the FC cooling protocol. Recently the 
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Fig. 6. t dependence of S at various T. t w = sec. H — 1 Oe. 

long-time decay of M^RAiit) has been also examined by 
Kenning et al. |27] for the same system with rapid FC 
cooling protocol and short wait time (t w = 7 - 110 sec). 
A ^-independent long-time decay overlaps with the in- 
dependent short-time decay. For the short-time decay, the 
corresponding S(t) exhibits a peak at t cr (w t w ). The long- 
time decay may be related to the initial state distribution 
developed during the FC cooling protocol. 

4.3 S(t) under the iJ-shift 

We have measured the t dependence of xzfc at T = 3.75 
K for various H, where t w = 1.0 X 10 4 and 3.0 x 10 4 sec. 
The measurements were carried out after the ZFC aging 
protocol: annealing of the system at T — 50 K and H = 
for 1.2 x 10 3 sec, quenching from 50 to 3.75 K at H = 0, 
and isothermal aging at 3.75 K for t w . The origin of t 
(t = 0) is the time just after H is turned on. The window 
time used in this measurement was 12 sec. Figures EJa) 
and (b) show the t dependence of S for t w = 1.0 x 10 4 
sec and 3.0 x 10 4 sec as H is varied as a parameter, where 
T = 3.75 K. The relaxation rate S(t) exhibits a peak at 
t = t cr , corresponding to a characteristic time scale of 
crossover from the isothermal aging state under H = 
to that under a finite H . Figure Eta) shows the plot of 
t cr at T = 3.75 K as a function of H for t w = 1.0 x 10 4 
and 3.0 x 10 4 sec. The value of t cr for t w — 1.0 x 10 4 
and 3.0 x 10 4 sec drastically decreases with increasing H . 
The value of t cr for t w = 3 x 10 4 sec is much larger than 
that for t w = 1 x 10 4 sec at low H (H < 50 Oe). How- 
(t)ever, they are almost identical irrespective of t w for high 
H (H > 100 Oe). We note that lnt cr is linearly dependent 
on H only for H < H a : H a w 50 Oe for t w = 3.0 x 10 4 sec 
and 100 Oe for t w = 1.0 x 10 4 sec. In Fig. Eft) we show 
the H dependence of S max . The peak value S max linearly 
decreases with increasing H irrespective of t w and tends 
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Fig. 7. t dependence of S for 5 < H < 300 Oe. T = 3.75 K. 
(a) t w = 1.0 x 10 4 sec. (b) t w = 3.0 x 10 4 sec. The time t = 
is a time when H is turned on. The ZFC aging protocol before 
the measurement is similar to that used in Fig. [5] 



to reduce to zero above 300 Oe. Similar behavior of t cr vs 
H has been reported by Zotev et al. [22] in Cu (1.5 at.% 
Mn). The break of the linear dependence of \xit cr on H 
occurs at H a . The value of H a decreases with increasing 
t w . Above H a the data of SM (= M TRM ~ M FC + M Z fc) 
vs H greatly deviates from zero, where Mtrm is the ther- 
mal remnant magnetization. This result suggests that the 
crossover from the quasi equilibrium to the noncquibrium 
regime occurs at H a , leading to the violation of FDT. 

According to Takayama |Hj, the H dependence of t cr 
under the 7J-shift aging process is governed by three lengths. 
LiT{t C r, H), Rx(t w ), and the crossover length Lh 1 given 

by 

L H /Lo « (H/Yh)- 1 / 6 , (7) 

where S — (d/2 — 9) and Th is the magnetic field corre- 
sponding to a wall stiffness T (a typical energy setting the 
scale of free energy barriers between conformations) . The 
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Fig. 8. H dependence of (a) t cr and (b) S ma x for t w — 1.0 x 10 4 
and 3.0 x 10 4 sec, which is obtained from Fig. [7] 



scaling relation is predicted to exist between the normal- 
ized lengths y = LT{t cri H = 0)/Lh and x = Rr(t w )/ 'Lh- 
y = x — chx 1+s with ch = 0.15 for x < 1. For i«l cor- 
responding to the case of low H and short t w , y — x, 
implying that t cr — t w . For x > 0.2 corresponding to 
large H and long t w , the curve y(x) deviates below the 
line y = x, indicating that t cr is shorter than t w . From 
this scaling relation, t cr can be obtained as 



t w [i - c H (H/r H )(t w /to) 5/z ^ T Y T) - 



(8) 



The following features of t cr vs H are derived from Eq. J3J) . 
(i) The time t cr is not simply proportional to t w . (ii) In 
the limit of H w 0, hit cr is linearly dependent on H: 



with 



ln(t cr /t w ) = -a H H, 



z(T)(c H /r H )(t w /t Q ) s /^ T l 



The slope an increases with increasing t w . (iii) The T de- 
pendence of t cr comes from the exponent z(T). The slope 
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olh increases with increasing T below T g mainly because 
of the factor (t w /t ) s/z ^ T) . 

We find that our results of Fig. IHJa) is consistent with 
the above predictions. The slope (an — 0.0505 ± 0.0003 
/Oe) for t w — 3.0 x 10 4 sec is larger than that (an — 
0.0262 ± 0.0005 /Oe) for t w = 1.0 x 10 4 sec. If we assume 
that 6 = 1.37, l/z{T) = bT/T g , T g = 3.92 K, b = 0.16, 
and t = t* = 5.29 x 10~ 6 sec [25|, then the value of T H 
can be estimated as Th — 2.15 kOe for t w = 3.0 x 10 4 sec 
and Th = 3.30 kOe for t„ = 1.0 x 10 4 sec. Our numerical 
calculation of ajj for 3 K < T < T g predicts that the slope 
an increases with increasing T below T g : oih — 0.036 
at T = 3.4 K and a H = 0.044 at T = 3.6 K for t w = 
3.0 x 10 4 sec. Although the measurement on an vs T has 
not been carried out in the present system, this prediction 
is consistent with the result reported by Zotev et al. [221 
in Cu (1.5 at.% Mn) as the increase of <xh vs T with 
increasing T for 0.7 < T/T g < 0.85. 

4.4 S(t) under the T-shift 

We have measured the t dependence of xzfc under the 
T-shift from the initial temperature T to the final tem- 
perature Tf (= 3.75 K = 0.957T S ), where = 3, 3.2 3.4, 
3.5, 3.6, and 3.9 K. The measurement was carried out after 
the ZFC aging protocol: quenching of the system from 50 
K to Ti, and isothermal aging at T — Ti and H = for t w 
(= 3.0 x 10 3 and 3 x 10 4 sec). The origin of time (t = 0) 
is a time just after T was shifted from T to Tf = 3.75 
K and subsequently H (— 5 Oe) was turned on. Figures 
EI a) and (b) show the t dependence of S at Tf = 3.75 
K at various initial temperature Ti. The temperature dif- 
ference is defined as AT = Tf — T L : the positive T-shift 
for AT > and negative T-shift for AT < 0. The re- 
laxation rate S(t) exhibits a peak at t = t cr irrespective 
of the sign of AT. The width of the peak in S(t) for the 
negative T-shift is much broader than that for the posi- 
tive T-shift. In Fig.^Jc) we show t cr as a function of AT 
for the positive T-shift for t w = 3.0 x 10 4 and 3.0 x 10 3 
sec. The decrease of t cr with increasing AT is observed at 
low AT for both t w = 3.0 x 10 4 and 3.0 x 10 3 sec. The 
value of t cr becomes independent of AT at large AT. For 
the negative T-shift, on the other hand, the value of t cr 
for AT = -0.15 K (T = 3.9 K) is larger than that for 
AT = OK (Ti = 3.75 K) for both t w = 3.0 x 10 4 and 
3.0 x 10 3 sec. Similar behaviors of t cr vs AT for the posi- 
tive T-shift have been reported by Granberg et al. [Cu (10 
at. %Mn)] [EH, Djurberg et al. [Cu (13.5 at. % Mn)] [H], 
and Jonsson et al. [Ag (11 at. % Mn)] [THj . 

The temperature chaos scenario postulates that the 
SG equilibrium configurations at different temperatures at 
Tf and T, are strongly correlated only up to the overlap- 
length Lat, beyond which these correlations decay to 
zero. From the droplet theory, the overlap length L^t 
is described by [T||2*] 

Lat/Lo « (T^lATl/Y^ 2 )- 1 ^, (9) 

where £ is the chaos exponent (£ = d s /2 — 8), d s is the 
fractal dimension of the surface of the droplet, and Tt 




L E ■ J . I E ■ " ■ ■ " 1 I I L 

■n- id* icr 

I (SEC) 



id- 



j a 1 1 ■ p 1 1 1 1 


1 J 1 1 | 1 1 1 1 | ■ 1 1 1 ^ 1 1 1 1 | 1 i. 1 1 , i 1 1 1 


\ 






i m = a ic io a sec 




— A 







E 0.1 OA 05 D.6 0? D.ft 



Fig. 9. t dependence of S at Tf = 3.75 K under the T shift 
from T to T f , where Ti = 3.0 - 3.9 K. (a) t w = 3.0 x 10 3 sec. 
(b) t w — 3.0 x 10 4 sec. The ZFC aging protocol is as follows: 
quenching of the system from 50 K to Ti, and isothermal aging 
at T = T and H = for t m . The measurement was started 
at t — 0, just after T was shifted from T to Tf = 3.75 K and 
subsequently H (= 5 Oe) was turned on. (c) t cr vs AT for 
the case of positive T-shift (t w = 3.0 x 10 3 and 3.0 x 10 4 sec), 
where Ti = Tf — AT (AT > 0) and Tf = 3.75 K. The solid 
lines are guides to the eye. 
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is the temperature corresponding to the wall stiffness Y. 
From this scenario a scaling relation is predicted to exist 
between the normalized lengths y = [LTf(t cr )/L^T] and 
x = [RTi{tw) / L at]'- y = x — ctx 1+ ^ with ct = 0.25 and 
C = 0.385 for x < 0.15 gUEE]. The cumulative aging 
corresponds to the relation y = x which is valid in the 
limit x w (the small \AT\ and short t w ). The large \AT\ 
and long t w corresponds to large x. For x > 0.05, the curve 
y deviates from the straight line y = x, corresponding to a 
rejuvenation due to the temperature chaos effect [211128) . 
The value of t cr /to can be described by 

t c r/t = (t w /t ) T ^ T ni-C T p(T f )\AT\(t w /t ) e/z{Ti Y Tf) , 

(10) 

with p(T f ) = Tf /2 /7%/ 2 . In the limit of \AT\ w 0, ln< cr is 
linear to T: 

\n(t cr /t w ) = -a T \AT\, 



with 



or 



= z(T f )c T p(T f )(t w /t f'^ 



In fact, the prediction that ax increases with increasing 
t w is experimentally confirmed from Fig. GUc) . The curve 
(Inter vs AT) is linearly dependent on AT (> 0) for AT < 
AT a with the slope a T - a T = 9.9 ± 0.4 and AT a = 0.3 K 
for t w = 3.0 x 10 4 sec and a T = 2.6 and AT a w 0.15 K for 
t w = 3.0 x 10 3 sec. If we assume that ( = 0.385, c T = 0.25, 
l/z(T f ) = bT f /T g , b = 0.16, and t = t* = 5.29 x 10~ 6 
sec |g] for t w = 3.0 x 10 4 sec, we havep(T/) = 1.62 ±0.07 
or T T = 1.12 ± 0.03 K for t w = 3.0 x 10 4 sec. 

Numerical calculations of Eq. (|10fl are carried out as a 
function of \AT\ for t w = t w (c&\) = 2.04 x 10 4 sec, T f = 
3.75 K, T t = T - AT, C = 0.385, t = r* = 5.29 x 10~ 6 
sec, and ct = 0.25. Note that we use t„,(cal) as t w in- 
stead of t w (exp) (= 3.0 x 10 4 sec) [^^(cal) = 0.68i tiI (exp)], 
which leads to the better agreement between our data and 
calculations at low \AT\. In Figs. HOf a) and (b) we show 
the results of t cr vs \AT\ for the positive and negative T- 
shifts, respectively, where p(Tf) is varied as a parameter 
from to 2.0. For comparison, our data of t cr vs AT with 
t w = 3.0 x 10 4 sec for the positive T-shift are also plotted in 
Fig. HUf al. For the positive T-shift, t cr decreases with in- 
creasing AT, independent of p(Tf) for < p(Tf) < 2. For 
the negative T-shift, t cr increases with increasing \AT\ for 
< p(Tf) < 0.7. It increases with increasing \AT\, show- 
ing a peak, and decreases with further increasing \AT\ 
for 0.7 < p(T f ) < 0.92 (see the inset of Fig. HCTb)). 
On the other hand, it decreases with increasing \AT\ for 
0.92 < p(Tf) < 2. This indicates that t cr decreases with 
increasing \AT\ for both the positive and negative T-shifts 
for 0.92 < p(T f ) < 2 (the symmetric T-chaos). We find 
that our data for both the positive and negative T-shifts 
agree well with the curve with p{Tf) w 0.8 in Fig. HUt a - ). 
corresponding to Yt = 1-80 K. This value of Yt is on the 
same order as that derived from the slope cut- 

For Ag (11 at. % Mn) |2T], t cr decreases with increasing 
AT for the positive T-shift. For the negative T-shift, how- 
ever, t cr shifts to the long-i side at \AT\ — 0.1 K and then 
t cr decreases with further increasing |^4T|. Similar behav- 
ior has been also observed in Cu (13.5 at. % Mn) |14| . 
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Fig. 10. Numerical calculations of t cr vs \ AT\ given 
for Tf = 3.75 K, T f = T + AT, and t w = t TO (cal) = 2.04 x 10 4 
sec, where T g = 3.92 K, b = 0.16, C = 0.385, c T = 0.25, 
t = t* = 5.29 x 10~ 6 sec, and < p(T f ) < 1. (a) The positive 
T-shift (AT > 0). For comparison, our data of t cr vs AT for 
t w — 3.0 x 10 4 sec are denoted by solid triangles. The value of t m 
(= t^(cal) = 2.04 x 10 4 sec) used in the calculation is different 
from that used in the experiment [t w (exp) = 3.0 x 10 4 sec], 
where U(cal) = 0.68t„(exp) (see Fig. 12b)). (b) The negative 
T-shift (AT < 0). Although the maximum of \AT\ is T g -T f = 
0.17 K in the present system, the value of t cr for the negative 
T-shift is independent of T g and is applicable to the other 
systems with different T 9 . 



where Tf = 58 K and t w = 3.0 x 10 3 sec. For the positive 
T-shift, t cr decreases with increasing AT (0 < AT < 1.21 
K). For the negative T-shift, t cr increases with increasing 
\AT\ for < \AT\ < 0.6 K, t cr becomes the longest at 
\AT\ = 0.6 K, and t cr decreases with further increasing 
\AT\ for 0.6 < \AT\ < 4.03 K. These behaviors can be 
well explained in terms of the scaling relations. 



4.5 Memory efFect in \' ar| d x" 

The memory effect in the SG system is defined as follows. 
When the system is cooled down to a low temperature 
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Fig. 11. T dependence of (a) A%' and (b) Ax"'- single mem- 
ory experiments (I and II). / = 0.05 Hz. h = 0.1 Oe. The 
measurement (I) (denoted as gradual) was carried out after 
the ZFC aging protocol (I): gradual decrease of T from 20 to 
3.75 K, isothermal aging at T\ = 3.75 K for 6.3 hours, and 
further gradual decrease from 3.75 to 1.9 K. The measurement 
(II) (denoted as quenched) was carried out after the ZFC ag- 
ing protocol (II): quenching of the system from 20 to 3.75 K 
at H — 0, isothermal aging at 3.75 K for 10 hours, and further 
quenching from 3.75 to 1.9 K. Both x' and x" were simultane- 
ously measured with increasing T (aging ZFC curves). Ax' and 
Ax" are the deviations of the aging ZFC curve from the refer- 
ence ZFC curve measured with increasing T after the standard 
ZFC aging protocol: quenching from 20 to 1.9 K at H = 0. 



below T g , a memory of the spin configurations which is 
imprinted in the specific cooling sequence, can be recalled 
when the system is re-heated at a constant heating rate 
|16II17| . In a single memory experiment, the memory is 
imprinted at T± for t w \ during the ZFC aging protocol. In 
a double memory experiment, the memory is imprinted at 
T\ for t w i and at T 2 (< T\) for t W 2 during the ZFC aging 
protocol. The dispersion and absorption thus recalled with 
increasing T are defined as x^ em (w,T) and x^ em (w,T), 
respectively. For comparison, the dispersion and absorp- 



tion as references [}(! re t{u,T) and x" e /( w :^)]> are also 
obtained with increasing T after the system is quenched 
from a high temperature above T g to the lowest temper- 
ature. Such AC susceptibility data are called as the ZFC 
reference susceptibilities, where no memory is imprinted. 
For clarity we define the difference between the aging 
ZFC and reference ZFC susceptibilities as Ax'(u>,T) = 
X'memi^, T)-x' ref {u, T) and as A X "(uj, T) = x'Lemi^, T)- 
x" e f(v,T). Figures UlT a) and (b) show the T dependence 
of Ax' and Ax" for the single memory experiment, where 
T x = 3.75 K, / = 0.05 Hz, and h = 0.1 Oe. Two kinds 
of measurements were carried out, depending on cooling 
rate during the ZFC aging protocol: (i) the gradual cool- 
ing (gradual decrease of T from 20 to 3.75 K at the cooling 
rate 5.6 x 10~ 4 K/sec, isothermal aging at T\ — 3.75 K 
for t w \ = 2.27 x 10 4 sec, and further gradual decrease 
from 3.75 to 1.9 K at the rate 2.6 x 10~ 4 K/sec), and 
(ii) the rapid cooling (quenching of the system from 20 
to 3.75 K at H = at the rate 0.13 K/sec, isothermal 
aging at T\ — 3.75 K for t w \ — 3.6 x 10 4 sec, and fur- 
ther quenching from 3.75 to 1.9 K at the rate 0.06 K/sec). 
Both x'mem and Xmem were simultaneously measured with 
increasing T (aging ZFC curves) from 1.9 to 8 K at the 
rate 1.7 x 10~ 4 K/sec. The reference ZFC curves {x' re f 
and x'ref) w ere measured with increasing T from 1.9 to 8 
K at the rate 1.7 x 10~ 4 K/sec after the standard ZFC 
aging protocol: quenching from 20 to 1.9 K at H = at 
the rate 0.15 K/sec. The results are as follows, (i) For the 
case of gradual cooling, Ax' exhibits negative local min- 
ima at 3.4 and 3.75 K, while for the case of rapid cooling, 
Ax' exhibits a negative local minimum at 3.75 K. Above 
4.3 K, Ax' is independent of the detail of the ZFC aging 
protocol, (ii) The absorption Ax" exhibits a negative lo- 
cal minimum at 3.75 K for both cases of the rapid and 
gradual cooling. Above 4.1 K, Ax" is independent of the 
detail of the ZFC aging protocol. In summary, the equi- 
libration at 3.75 K gives rise to a dip in Ax' and Ax" , 
suggesting that the memory of spin configurations which 
are imprinted during the cooling process, is recalled dur- 
ing the heating. The dip of Ax' and Ax" at 3.75 K for 
the rapid cooling is much narrower than that for the grad- 
ual cooling. The reason is that in the case of the gradual 
cooling the spin configurations at T not equal to 3.75 K 
are also imprinted during the cooling process. Similar but 
more pronounced single memory effect has been reported 
by Jonsson et al. for Ag (11 at. % Mn) [TH] , 

Figures IT2T a) and (b) show the T dependence of Ax' 
and Ax" for the double memory experiment, where T\ = 
3.75 K, T 2 = 3.0 K, / = 0.05 Hz, and h = 0.1 Oe. The 
measurement was carried out after the ZFC aging proto- 
col (only in the case of rapid cooling): quenching of the 
system from 20 to 3.75 K at H = at the rate 0.13 K/sec, 
isothermal aging at T% = 3.75 K for t w \ = 3.6 x 10 4 sec, 
quenching from 3.75 to 3.0 K at the rate 0.025 K/sec, 
isothermal aging at T2 = 3.0 K for t W 2 — 3.6 x 10 4 sec, 
and quenching from 3.0 to 1.9 K at the rate 0.04 K/sec. 
Both x'mem and Xmem were simultaneously measured with 
increasing T (aging ZFC curves). Ax' and Ax" are the de- 
viations of the aging ZFC curves from the reference ZFC 



M. Suzuki, I. S. Suzuki: Aging, rejuvenation, and memory effects in Cuo.sCoo.sC^-FeCLj GBIC 



11 



0.02 



CU 0.5 C °0.5 CI 2- FeCI 3 GBIC 



CU 0.5 CO 0, CI 2- FeCI 3 GBIC 



> 
E 

5 



-0.02 



-0.04 



-0.06 



-0.08 



> 

1 -0.005 

a). 



-0.01 



-0.015 





(a); 






• 

• • 


• 

• 






_ • • 




... 






• 


- 

• 


• 


T - ^ 7R K 
i , — o. / \j r\ 


- 

- 


T 2 = 3.0 K 


- 


• 




1 1 M 1 M 1 1 1 1 1 M 1 1 1 


1 1 1 1 1 1 | 1 1 1 1 1 1 1 1 

. • (b) : 

• 


- 

• 

• 

_ • 
• 


• 

• • • " 










-• • 






• 


• 

• 




•• • 








• 




• 









2.5 



3.5 4 
T(K) 



4.5 5 5.5 



Fig. 12. T dependence of (a) Ax' and (b) Ax"' double mem- 
ory experiment. / = 0.05 Hz. h = 0.1 Oe. H — 0. The measure- 
ment was carried out after the ZFC aging protocol: quenching 
of the system from 20 to 3.75 K at H = 0, isothermal aging at 
3.75 K for 10 hours, quenching from 3.75 to 3.0 K, isothermal 
aging at 3.0 K for 10 hours, and quenching from 3.0 to 1.9 K. 
Both x an d x" were simultaneously measured with increasing 
T (aging ZFC curves) . Ax' and Ax" are the deviations of the 
aging ZFC curve from the reference ZFC curve measured with 
increasing T after the standard ZFC aging protocol: quenching 
from 20 to 1.9 K at H = 0. 



curves. The results are as follows. Both Ax' and Ax" show 
a small dip at T 2 = 3.0 K and a large dip at 4.2 K. There 
is no dip at T\ = 3.75 K. The magnitudes of Ax' and 
Ax" for the double memory effect are smaller than those 
for the single memory effect. It seems that the spin con- 
figurations which are imprinted at T\ = 3.75 K during the 
cooling process may be partially reinitialized by the spin 
configurations imprinted at Ti — 3.0 K. Similar (but more 
pronounced) double memory effect has been reported by 
Jonsson et al. for Ag (11 at. % Mn) ^HJ. The condition 
for the appearance of the dips at T\ and T 2 (= T\ — AT) 
is that the overlap distance Lat given by Eq.© with 
T = T2 is larger than the size Rrtitwi)- The spin con- 
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Fig. 13. Relaxation of x'( w >*) and x"(u,t) at T 2 = 3.55 K 
during a temperature cycle between Ti = 3.75 K and T2 = 3.55 
K. The change of T with t is also shown, t — is a time just 
after the standard ZFC aging protocol: quenching of the system 
from 50 to 3.75 K at H = 0. The data of x"(w,t) are the same 
as those presented in Ref. I2(i|. 



figurations imprinted at T% is partially reinitialized at T2, 
when AT is larger than the threshold temperature (AT) t 
given by 

(AT) t = (T^ /2 /Tl /2 )(t wl /t )-^ z ^\ (11) 



When f. 



3.6 x 10 4 sec, t Q = t* 



5.29 x 10" 6 sec, and 

„3/2 



C = 0.385, {AT) t is estimated as {AT) t = 0.147^". Using 
the value of T T (w 1.12 - 1.8 K) obtained in Sec. El the 
value of (AT) t is estimated as (AT) t = 0.16 - 0.33 K, 
which is smaller than AT = 0.75 K. 



4.6 Rejuvenation effect in x' ar| d x" under the T- and 
/T-shifts 

The rejuvenation and memory effects in x' an( l x" un_ 
der the T-shift are observed in our system. Figures H~3T a) 
and (b) show the t dependence of x'(cj,i) and x"(w,t) at 
Ti = 3.55 K under the T-shift between Ti = 3.75 K and 
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T 2 = 3.55 K, where / = 0.05 Hz and h = 0.1 Oe. Here 
the data at T\ — 3.75 K are not shown (see the data of 
X" at T\ = 3.75 K in the previous paper [2H|)- First our 
system was quenched from 10 K to T\ = 3.75 K at H = 0. 
The origin of t (t = 0) is a time when T becomes T\ . The 
relaxation of x' and x" was measured as a function of t 
during a period t w i 8.2 x 10 3 sec). The temperature 
was then changed to T 2 (the negative T-shift). The relax- 
ation of x' and x" was measured as a function of t for a 
period t W 2 (~ 8.2 x 10 3 sec) at T 2 . The system was again 
heated back to T\ (the positive T-shift). These processes 
were repeated subsequently. Just after every negative T- 
shift, both x' and x" do not lie on the reference curves of 
X' and x" at T 2 obtained when the system is quenched to 
T 2 directly from 10 K at t = 0. The values of x' and x" 
are larger than the reference curves, indicating the partial 
reinitialization (rejuvenation) in x' and x" after the neg- 
ative T-shift. Just after every positive T-shift, however, 
both x' and x" ne on to the reference curves of x' and x" 
at Ti obtained when the system is quenched to T\ directly 
from 10 K at t = 0, indicating the memory effect in x' 
and x" after the positive T-shift. Note that the reference 
curve coincides with a curve where the lowest points for 
each relaxation in x' and x" are connected as a function of 
t. The strong rejuvenation effect for the negative T-shift 
is also predicted from numerical study by Takayama and 
Hukushima |Zj using the MC simulation on the 3D Ising 
EA SG model. Here we note that the threshold temper- 
ature differences under the positive and negative T-shift 
between Ti and T 2 , {AT) + and (AT)., are given by 

(AT) + = (r^7T 2 1/2 )(W*or CMT2 \ 

(AT). = {T^ 2 /T 1 1 ,2 ){t wl /t )-^ z ^\ (12) 

respectively. When Ti = 3.75 K, T 2 = 3.55 K, t wX = t w2 = 
8.2 x 10 3 sec, to = t* = 5.29 x 10~ 6 sec, and C = 0.385, 
we have (AT) + = 0.148T, 3 / 2 and (AT). = 0.163T, 3 / 2 : 
(AT)+ = 0.18 K and (AT). = 0.19 K for T T = 1.12 K. 
Thus the values of (AT) + and (AT), are comparable to 
the difference T x - T 2 = 0.20 K. 

The rejuvenation effect in x' and x" under the if-shift 
is also observed in our system. Figures IllT a) and (b) show 
the t dependence of x' and x" at T = 3.75 K under the 
F-shift, where / = 0.05 Hz and h = 0.1 Oe. After the ZFC 
aging protocol (quenching from a high temperature above 
T g to a temperature T (= 3.75 K) ), x' and x" at H = 
were measured for a period t w i (w 1.3 x 10 4 sec), where 
the origin of t (t = 0) is a time when T becomes 3.75 K. 
The field is changed from to H (= 50 Oe) at t — t w %, 
After this 7J-shift, x' aim x" were measured for a period 
t W 2 (k 1.3 X 10 4 sec). Subsequently, the field was turned 
off from H to and the measurements were carried out at 
H = for a period t w %. This process was repeated. We find 
that both x' and x" undergo drastic jumps under the H- 
shift (from to 50 Oe) to values higher than the reference 
curves at H = 50 Oe. They also undergo drastic jump 
under the i?-shift (from 50 to Oe) to values higher than 
the reference curves at H = 0. These results suggest that 
a partial reinitialization (rejuvenation effect) of x' and x" 
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Fig. 14. Relaxation of x (to, t) and x"^! t) during a magnetic- 
field cycle between H — and 50 Oe. T = 3.75 K. t = is a 
time just after the standard ZFC aging protocol: quenching of 
the system from 50 to 3.75 K at H = 0. The data of x"( w 7*) 
are the same as those presented in Ref. |2(i| . 



occurs for the iTshifts (0 to 50 Oe and 50 to Oe). Note 
that the point at T = 3.75 K and H as 50 Oe is located on 
the de Almeida-Thouless (AT) [23 line in the (H,T) phase 
diagram [SB]. On this AT line, the T-derivative d<5/dT 
shows a local minimum, where S — xfc ~ Xzfc- 

When the overlap distance Lh given by Eq. is smaller 
than the size Rt^w), the spin configuration imprinted at 
H = is partially reinitialized at H (> H t ), where the 
threshold magnetic field is defined as 



Hi 



r H (t w /h)- s '< T \ 



(13) 



When t w = 3.6 x 10 4 sec, t = t* = 5.29 x 10~ 6 sec, 
S = 1.37, and T = 3.75 K, H t is estimated as H t = 8.67 x 
10~ 3 Xh . Using the value of Th (=2-3 kOe) obtained in 
Sec. 14.31 the value of H t is estimated as 17 - 25 Oe, which 
is lower than 50 Oe. 
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5 Discussion and conclusion 

The non-equilibrium nature of the spin dynamics in 3D 
Ising SG Cuo.sCoo.sC^-FeCla GBIC has been studied from 
the t dependence of xzfc, x' > and x" after specific ZFC 
protocols including the T- shift (AT) and if-shifts. The 
relaxation rate S(t) shows a peak at t cr (ss t w ), corre- 
sponding to a crossover from quasi-equilibrium dynamics 
to non-equilibrium dynamics. The value of t cr strongly 
depends on the wait time t w , T, H, and AT. The rejuve- 
nation effects are observed in x' and x" under the negative 
T-shift and both the positive and negative if -shifts. The 
spin configurations imprinted under the ZFC aging proto- 
cols are recalled on heating the system. 

We find that the observed change of t cr under the 
T-and ii-shifts is well explained in terms of the scaling 
relations where the overlap lengths Lat and Lu play 
a significant role. Under the T-shift from T = T to Tf 
(= T + AT), the size of domains are unaffected for suf- 
ficiently small | |, where Lat is larger than RTi{t w ). 
Then the relaxation rate S(t) shows a peak at t cr where 
LiTf (t) is equal to (t w ). In contrast, the size of domains 
are affected for sufficiently large \AT\. The overlap length 
Lat becomes lower than RTi(t w )- Then S(t) has a peak 
at t cr , where LT f {t) is equal to Lat- Under the iJ-shift, 
the size of domains are unaffected during the H shift for 
sufficiently small H, where the overlap length Lh is larger 
than RT(t w ). The value of t cr is dependent on t w . In con- 
trast, the size of domains are affected during the H shift 
for sufficiently large H, where Lu becomes lower than 
Rritw)- Then S(t) has a peak at t cr , where Lx(t) ~ Lh- 

We discuss the scaling relation of the T dependent re- 
laxation rate S(T, t). As pointed out in Sec. l4.1l YzFr'fr. t) 
below T g may be described by a scaling function 

XzFc{T,t) = G(x), (14) 

with x = LT{t) / Rr{t w ). From the definition, S(T,t) is 
derived as 

S(T,t) = d X zFc(T,t)/dliit = (l/z(T))H(x), (15) 

with H(x) = x(dG(x)/dx). It follows that S(t) is de- 
scribed by a scaling function H(x) except for the factor 
l/z(T). When the scaling function H(x) has a peak at 
x = a (constant), then S(t) has a peak at t = t cr , where 
x = LT(t) / Rr{t w ) = a. Then t cr is simply described as 

\n{t cr lt w ) = (T g /bT)\na. (16) 

When a is larger than 1, hi(t cr /t w ) increases with decreas- 
ing T. In fact, the least squares fit of the data of ln(t cr /t w ) 
vs t with t w = 2.0 x 10 3 sec in Fig. OJa) (H = 1 Oe) to 
Eq.OH) yields to (T g /b) In a = 37.1 ± 1.9. Since T g = 3.92 
K and b = 0.16, we have a = 4.5 ± 0.4. It is predicted 
from Ea. l|15J) that S max increases with increasing T for 
T < T g since S max is linearly dependent on the factor 
l/z(T) (= bT/T g ). Experimentally, as shown in Fig.^b) 
with t w — 2.0 x 10 3 sec, the peak height S max strongly de- 
pends on T. The peak height S max exhibits a broad peak 
around 3.6 - 3.8 K, just below T g . The linear increase of 



Smax with increasing T below T g is considered to be due 
to the factor 1/ z(T). The decrease of S max with increasing 
T above T g , however, cannot be explained in terms of the 
above model because the scaling relation is valid only for 
T < T g . In Sec. Owe show that S(t)/S max at T = 3.75 
K and H = 5 Oe obeys a t / '^-scaling law for long t w 
(5.0 x 10 3 < t < 1.5 x 10 4 sec) and < t < 6.0 x 10 4 sec 
(see Fig. EIc)): S(t)/S max is well described by a scaling 
function F(t/t w ) which has a peak at t cr /t w sa 0.68. Using 
Ea. (|16(l . the value of a is estimated as a = 1.11, indicating 
that the constant a at H = 5 Oe is different from that at 
H = 1 Oe. 

Finally we consider the cause of the complicated be- 
havior of S(t) at T = 3.75 K for < t w < 750 sec. One 
of the cause is the way how T approaches 3.75 K during 
the ZFC aging protocol. The temperature drops rapidly 
to « 3.50 K and slowly approaches 3.75 ± 0.01 K from the 
below (usually) within 230 sec. The experiment of S(t) 
under the T shift (see Sec. 14.4(1 suggests that the initial 
undercool is not so important because of the temperature 
difference larger than the threshould temperature differ- 
ence (AT) t [231 - However, the subsequent approach of T 
to 3.75 K for a wait time t W Q (< 230 sec) may play a signif- 
icant role in the aging behavior after t =0 at T m — 3.75 K. 
The threshold temperature difference (AT) t is estimated 
as 



(AT) t = (7^ /2 /[(T Tn -(^T) t ) 1 /2])( tto0 /t )-C/*(T m -(2 1 D t ) ) 

from the condition that Tat = Rr-Aritwo)- When to = 
t* = 5.29 x 10~ 6 sec, C = 0.385, and t w0 = 230 sec, 
(AT)t is estimated as (AT) t = 0.24 K for T T = 1.12 K 
and 0.56 K for Tt = 1.80 K. This suggests that the spin 
configuration at T m = 3.75 K after t = is affected by 
that imprinted at T m — AT at the wait time t W Q. Note that 
the value of t cr at T m — AT is larger than that at T m for 
the same t w (see Fig. CHa)). The appearance of two peaks 
in S vs t at < t w < 200 sec may be associated with two 
domains generated for the wait time t w o at T m — AT and 

fOr tyj at Tyyi . 

In conclusion, we have undertaken an extensive study 
on the t dependence of the relaxation rate S(t) mainly 
below T g under the various conditions including the T- 
and i/-shifts. The t dependence of S is well explained in 
terms of the scaling relation. The peak of S(t) occurs when 
the mean SG domain size Lx(t) coincides with RT(t w ), 
where t w is the wait time. Our results indicate that the 
aging, memory and rejuvenation phenomena observed in 
our system are very similar to those in conventional spin 
glass systems. 
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